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single simulation. 
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I. INTRODUCTION 

With the rapid development of computer processors, 
numerical simulations using the Monte Carlo method 
have become a well established tool for the study of pro- 
teins, polymers and spin-glass models. When a Monte 
Carlo simulation using Metropolis importance sampling 
is carried out at a fixed temperature, the quality of 
the data depends on the distance from criticality, and 
therefore, multiple computer runs should be performed 
for each temperature. In order to speed up the simulation 
several methods have been suggested to overcome prob- 
lems such as slow dynamics which makes the Metropolis 
algorithm inefficient, and to study systems with a rough 
energy landscape with multiple local minima in free en- 
ergy. Important examples are the cluster-flip Swendsen- 
Wang algorithm ;2j, which has been used to reduce crit- 
ical slowing down at second order phase transitions, the 
multicanonical method [j| which was introduced to over- 
come the barrier between coexisting phases at first order 
phase transitions, the broad histogram method j4<] which 
directly calculates the density of states with only one 
computer run, and the flat histogram method 5] which 
generates a flat histogram in energy space similar to mul- 
ticanonical simulations. The drawback of the flat his- 
togram method is the slow diffusion of the random walk 
which is the same as in the multicanonical method. Nev- 
ertheless, no method is more efficient than that recently 
proposed by Wang and Landau [E El which allows 
one to get around these difficulties even for large systems. 

In this work we present results of simulations of three- 
dimensional lattice polymer chains ^(j using Wang- 
Landau sampling, and calculate thermodynamic proper- 
ties of the system. The method is described in Sec. II. 
Sec. Ill provides a brief background on the algorithm. 



The results are given in Sec. IV and we summarize and 
conclude in Sec. V. 



II. THE MODEL 

We consider a homopolymer consisting of N monomers 
which may assume any self avoiding walk (SAW) config- 
uration on a three-dimensional lattice. Assuming that 
the polymer is in a bad solvent, there is an effective 
monomer-monomer attraction in addition to the self- 
avoidance constraint representing excluded volume. For 
every pair of nonbonded nearest-neighbor monomers the 
energy of the polymer is reduced by e (See Fig. QJ. 




FIG. 1: Interaction Self Avoiding Walk. 

Therefore, the interaction between two nonbonded 
monomers i and j is given by 



* Corresponding author: 
cunha@posgrad. ufg. br 
t caparicaOif . ufg. br 



where 6 is the lattice constant. 
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The Hamiltonian for the model can be written as 

n = -eJ2 ( 2 ) 



<i,3> 

where a = 1 (0) if the site i is occupied (vacant), and 
the sum is over nearest- neighbor pairs (The sum 

is understood to exclude pairs of bonded segments along 
the chain). We are interested in the equilibrium proper- 
ties at temperature T. (In the following, temperature is 
measured in units of e/fce, where fcs is the Boltzman's 
constant.) 

In order to generate a Markov process for sampling 
configurations we have first to define a protocol of move- 
ments. We used the so-called reptation or "slither- 
ing snake" move which consists of randomly adding 
a monomer to one end of the chain and removing a 
monomer from the other end, maintining the size of the 
polymer constant. We define one Monte Carlo step as N 
attempted moves. At certain moments, one end of the 
chain may not be able to move, but successive motions 
of the other end release it from the trap. 



III. THE METHOD 

The Wang-Landau algorithm calculates the density 
of states g{E) by carrying out a random walk in 
energy space with an acceptance probability propor- 
tional to l/g(E) instead of the usual Boltzmann weight 
e -E/k B T usec j m conventional Monte Carlo simula- 
tion. The probability of energy E changes therefore from 
g{E')e~ E l kBT in the Metropolis scheme to an almost con- 
stant probability in the Wang-Landau method. As a re- 
sult the probability minimum which appears at first order 
phase transitions, for example, practically vanishes. The 
simulation is performed such that if E\ is the energy of 
the current configuration and Ei is the energy of a possi- 
ble new configuration, the acceptance probability is given 
by 



p{E x 



E2) = rain 



g{Ei) 
g{E 2 ) 
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(3) 



For each accepted configuration we accumulate an energy 
histogram H(E). Since the density of states is not known 
a priori, Wang and Landau proposed to set g(E) = 1 
initially, for all energy levels. 

To study large systems the energy space is divided into 
bins and the random walk is carried out independently 
in each segment. As the random walk in energy space 
is performed, whenever a move to a configuration with 
energy E is accepted the density of states is updated 
by multiplying it by a "modification factor" / > 1 that 
accelerates the diffusion of the random walk, and an unit 
is added to the histogram H(E), i.e. 

g(E) -> g(E) ■ /; and H(E) — > H(E)+1. (4) 



FIG. 2: Density of states for a chain of N = 50 monomers. 
The ground state of this chain is E g — —58. 



The initial choice of / is / = e = 2.71828.... The den- 
sity of states is multiplied by / until the accumulated 
histogram H(E) is flat. We then reduce / by setting 



(5) 



and resetting H(E) = for all energy values. 

This process is repeated; the simulation converges to 
the true value of g{E) when / is approximately 1. In our 
simulations the criterion of flatness was taken as each 
value of the histogram reaching at least 80% of the mean 
value (H(E)). The histograms are generally checked af- 
ter each 1000 Monte Carlo steps. 



IV. RESULTS 

In this section we present results obtained for chains of 
lattice homopolymers using the Wang-Landau algorithm 
and compare them with those obtained using Metropolis 
algorithm importance sampling. In Fig. [5] we show the 
density of states for a chain of 50 monomers. 

The ground state is achieved by shaping a configura- 
tion rolled up like a snail (see Fig. |3J). The energy of 
this configuration is obtained by counting the number 
of nonbonded nearest-neighbor pairs. In our simulations 
the polymer was always set initially in this configuration 
and, using the algorithm described above, the density of 
states was estimated as a result of a random walk in the 
energy space. 

Knowing the density of states, one can calculate any 
thermodynamic property A of the system through the 
canonical average 



(A) T = 



J2 E (A) E g(E)e- E / k - T 



(0) 



where (A) e is the microcanonical average of observable 
A obtained during the simulation, g{E) is the density 



3 




Metropolis 





FIG. 3: A ground state configuration for N = 50. In this case 



the energy is E g 



-58. 



of states for each energy level E, and Z is the partition 
function 



Z 



-E/k B T 



(7) 



In Monte Carlo simulations using Metropolis sampling 
we fix the temperature and determine the thermody- 
namic average {A)t- Therefore a new simulation is 
needed for each value of temperature. In Wang-Landau 
sampling we estimate the density of states and then cal- 
culate any thermodynamic property of interest by means 
of a simple algebraic operation. The computational time 
in our simulations was about ten times smaller using the 
Wang-Landau method as compared with the Metropolis 
algorithm. 

In Fig. 0|we show the result for the internal energy via 
Wang-Landau sampling, evaluated using the relation 



-E/k B T 



(S) 



and compare it with that obtained via the Metropolis 
algorithm. 

The specific heat can be determined from the fluctua- 
tions in the internal energy 



C V (T) 



(E 2 ) T ~ (E)\ 
k B T 2 



(9) 



The temperature dependence of the specific heat is 
shown in Fig. |SJ Much more CPU time would be needed 
to obtain significantly better results using the Metropo- 
lis algorithm. One can see clearly from these results that 
the Wang-Landau method yields a better description in 
the low-temperature regime. 

The temperature dependence of the mean square end- 
to-end distance, given by 

(R 2 ) = ([{x N - x Q ) 2 + {y N - y ) 2 + (z N - z ) 2 }) (10) 



FIG. 4: Temperature dependence of the internal energy for a 
chain of N = 50 monomers. 
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FIG. 5: Specific heat vs. temperature for a chain of N = 50 
monomers. 
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FIG. 6: Mean square end-to-end distance vs. temperature for 
a chain of iV = 50 monomers. 



is shown in Fig. H3 

Differentiating expression with regard to tempera- 
ture using (R 2 ) as the observable A, we obtain the fluc- 
tuation of the mean square end-to-end distance as 



d{R 2 ) T {ER 2 ) T -{E) T {R 2 ) rj 



dT 



2^2 



(11) 



4 



Metropolis 




Wang-Landau * k 




FIG. 7: Derivative of the mean square end-to-end distance vs. 
temperature. The curve with a peak resembles a second-order 
phase transition. 
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cannot be calculated directly by conventional Monte 
Carlo simulation. It can be estimated by integration over 
other thermodynamic quantities, such as specific heat, 
but these calculations are somewhat unreliable since the 
specific heat itself is not easy to estimate accurately 
With an accurate density of states estimated by the 
Wang-Landau method, the entropy can be calculated 
simply using 



S(T) = 



E(T) - F(T) 
f ' 



(13) 



In Fig. [!|]we show our result for the temperature depen- 
dence of the entropy for a homopolymer chain of N — 50 
monomers. 



FIG. 9: Entropy vs. temperature for a chain of iV = 50 
monomers. 



FIG. 8: Free energy vs. temperature for a chain of N = 50 
monomers. 

The temperature dependence of the derivative of the 
mean square end-to-end distance is shown in Fig. [7] The 
coil-globule phase transition is characterized by a peak 
in the derivative of (R 2 ) similar to the specific heat. We 
believe that the low-temperature peak corresponds to ex- 
citations of the surface of the typlically compact, folded 
configuration. This effect is more pronounced in short 
chains. 

An advantage of this algorithm is that we can readily 
calculate the free energy and entropy, quantities that are 
not directly accessible in conventional Monte Carlo simu- 
lations. In terms of the density of states, the free energy 
can be calculated from the partition function 

F(T) = -k B T\n(Z) = -k B T\n g(E) e - E ^ T ^ . 

(12) 

The temperature dependence of the free energy ob- 
tained from our simulation is shown in Fig. [S] for a ho- 
mopolymer with N = 50 monomers. 

The entropy is a key thermodynamic quantity that 



V. CONCLUSIONS 

We study the thermodynamic behavior of a three- 
dimensional homopolymer chain in lattice using the 
Wang-Landau sampling. We show that the density of 
states directly obtained by this algorithm enables one 
to calculate thermodynamic properties even for large 
systems with only one computer run. We also note 
from the temperature dependence of the specific heat 
and the mean square end-to-end distance that the low- 
temperature region is better explored by this method 
than by standard Monte Carlo simulation. Finally, ther- 
modynamic quantities such as the free energy and the 
entropy, which are not easily accessible by conventional 
methods, are directly estimated from the density of 
states. 
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